Effect of resonances on the transport properties of two-dimensional disordered systems 
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We study both analytically and numerically how the electronic structure and the transport 
properties of a two-dimensional disordered system are modified in the presence of resonances. The 
energy dependence of the density of states and the localization length at different resonance energies 
and strengths of coupling between resonances and random states are determined. The results show, 
that at energy equals to the resonance energy there is an enhancement in the density of states. 
In contrast, the localization length remains unaffected from the presence of the resonances and is 
similar to the one of the standard Anderson model. Finally, we calculate the diffusion constant as a 
function of energy and we reveal interesting analogies with experimental results on light scattering 
in the presence of Mie resonances. 
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The study of waves propagating in disordered lattices, 
during the past years has attracted much attention and 
many interesting phenomena have been, more or less, well 
understood. However, until recently the interest of the 
physics community was primarily focusing on quantum, 
i.e. electronic waves. It was the observation of the co- 
herent backscattering effect in classical wave systems , 
the analogous effect to weak localization in the electronic 
case, which triggered a burst of interest in further stud- 
ies of scattering disordered classical wave systems. The 
question of localization of classical waves has attracted 
attention for two reasons. First, the properties of classi- 
cal waves such as light waves, microwaves, and acoustic 
waves in random media are of fundamental interest for 
their own sake. Second, classical waves can serve as a 
model system for testing the theory of Anderson localiza- 
tion of electrons experimentally in a clean way, without 
the complication of strong inelastic scattering and other 
effects of electron- electron and electron-phonon interac- 
tion. Existing theories predict the localization of classical 
waves under certain circumstances ||l[ leading thus to the 
conclusion that the analogy between quantum and clas- 
sical waves works reasonably well. However, there is no 
conclusive experimental evidence, and thus the full cor- 
respondence is not yet established beyond any doubt (a 
thorough discussion may be found in Refs. 

The present paper has been motivated by an experi- 
mental work, performed several years ago Q on the scat- 
tering of light by a disordered medium in the presence of 
Mie resonances. A low concentration of Mie resonances 
leads to a strong change in the transport properties of the 
system, reflected in a strong reduction of the transport 
velocity, or equivalently, of the diffusion constant D{E), 
near the resonance energies. With increasing concentra- 
tion of resonant scatterers, the dip in D{E) widens and 
becomes less deep. Qualitatively speaking, the transport 
velocity is reduced because a lot of energy is temporar- 



ily stored inside the resonance or equivalently the wave 
spends a lot of time (dwell time) inside the Mie reso- 
nances. 

To understand this surprising result, different theoreti- 
cal approaches based on different considerations were de- 
veloped [Q, However, there are controversies as 
to whether this can be applied to the quantum case, i.e., 
electronic waves. In fact, at low concentration of Mie 
resonances Q , in the spirit of the coherent-potential ap- 
proximation (GPA) 0], it was shown that, unlike elec- 
tronic systems, the diffusion constant D{E) for classical 
waves decreases sharply close to the resonance energies. 
This result was viewed as due to the different Ward iden- 
tity caused by an energy-dependent scattering potential 
in the classical case. An extension of the CPA to the 
strong concentration limit the so called coated CPA, 
shows that the effect of the Mie resonances decreases at 
this regime. On the other hand, in a recent work ||^, it 
was shown that an extension of the Random Matrix The- 
ory (RMT) [0 , capable of accounting for the presence of 
resonance scattering, is able to explain the obtained ex- 
perimental result and generalize it to both quantum and 
classical waves, although the corresponding Ward identi- 
ties are different 

In the present paper, we extend our previous anal- 
ysis on the effect of resonances on electrons in a one- 
dimensional disordered medium to higher dimensions. 
The quantities that will monopolize our interest, will be 
the density of states (DOS) p{E) and the localization 
length loo{E). In the one-dimensional case at low 
concentration and weak coupling both these quantities 
are affected strongly by the presence of resonances. The 
DOS exhibits a Lorentzian peak at the resonance energy 
and the localization length is drastically decreased except 
at the resonance energy which remains essentially the 
same with the one corresponding to the standard Ander- 
son model. We will show that in higher dimensions, the 
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density of states exhibits the same resonant enhancement 
while the locahzation length remains unafFected in con- 
trast to the one-dimensional case. Finally, by making use 
of the well known results for quasi-one-dimensional sys- 
tems we will compute the diffusion constant D{E) and 
identify similarities with the corresponding problem of 
light scattering in the presence of Mie resonances. 

The mathematical model we consider is the tight- 
binding Anderson Hamiltonian on a two-dimensional lat- 
tice, 



H 



ei\l >< l\+^tu,\l X I' 



(1) 



where (/) denote the sites of a two-dimensional lattice 
L X M . The local site energies e/ are taken indepen- 
dently at random within the interval [—Wl2,W/2]. At 
some randomly chosen sites k, (m in total), the energies 
ek are taken to be fixed i.e. — E^- These energies cor- 
respond to the resonance energies. To simplify the prob- 
lem, we assume that all resonances have the same energy 
Er- The hopping matrix elements tui are restricted to 
nearest neighbors. When the coupling is between ran- 
dom states {N in total), tw are constant taken to be 
the unit of energy while the coupling between a resonant 
state and another state tw has a different constant value 

V. 
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FIG. 1. The density of states p(E) for the 2d-Andersfii model in 
the presence of resonance states with energy Er = and various 
values of the coupling constant ii corresponding to different sym- 
bols. The thick dashed line corresponds to the Anderson model 
without resonances. In all cases the disorder strength is W = i. 

To evaluate the density of states p{E) of this system, 
we start by writing the original [N +m) x +m) Hamil- 
tonian matrix H in the following form 



H 



Ho V 

Er X /„ 



(2) 



This form can be obtained easily by separating the res- 
onant states from random states. V"^ denotes the trans- 
pose of V and Hq is a N- dimensional matrix that 



describes the part of the Hamiltonian H without reso- 
nances. Here Im is the m-dimensional unit matrix, and 
Er is the resonance energy of each of the m resonant 
scatterers. The rectangular matrix V couples the m res- 
onances to Hq. All the matrix elements of each column 
of V are zero except four with value v. 
The Green's function corresponding to the Hamiltonian 
(0), can be written as 



G = 



Go 

Gy 

E^ 



Gy 
Gr 



ij] - Hq 



-V 

{E - Er) X /„, + ir] 



(3) 



To determine the different matrices in G, we use the pre- 
vious relation to obtain the following set of equations: 



{E-Ho + iT])Gv - VGr = , 
-V^Gv + [{E - Er)I + iri]Gr ^ 1. 

Solving (Q) with respect to Gr we get the expression 
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G, = 



{E - Er)I + lll-VTGaV 



(4) 



(5) 



where Gq = -^zrw^ppi:^- In the limit of weak coupling 
u <C 1 and low concentration of resonances one can make 
the approximation Gq ~ Gq. Then the density of states 
is given by 



(6) 



which, in the limit of u ^ 1 and low concentration of 
resonances, can be written as 



p{E)^po{E) 



r/27r 



{E 



(7) 



where po{E) ~ — ^TrGg. To a good approximation, 
Po{E) can be identified as the density of states of the 
standard two-dimensional Anderson model without res- 
onances. The effect of the resonances is given by the 
second term on the right hand side of Eq. (|^). It ex- 
hibits a resonance enhancement near the resonance en- 
ergy Er described by a Lorentzian peak. The width, 
r « 2TrpQv'^, is given by the decay width of individual 
resonances. Based on physical grounds, we expect that 
in the limit of weak coupling (w ^ 1) there will be a 
strong degeneracy due to the resonance states. This de- 
generacy is responsible for the resonance enhancement of 
the DOS at E — Er- Gradually, as the coupling with 
the other states is increased, the degeneracy and conse- 
quently the peak in p{Er) decrease. These expectations 
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are in full agreement with Eq. (^. The same result has 
been found in the extended RMT, Ref. with the use 
of the super-symmetry formalism. Here however Eq. (|^) 
appears natural under the frame of the tight-binding pic- 
ture. 
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concentration / of resonances in the disordered sample. 
In the weak coupling limit, we found that by increas- 
ing the concentration / the Lorentzian peak at the DOS 
increases in high (see Fig. 2) in agreement with the ana- 
lytical formula (^. Our numerical calculations, for small 
resonance coupling v ~ 0.2, are presented in Fig. 2. for 
two particular values of concentration. We remind that 
exactly the same behavior of p{E) was observed for the 
one-dimensional case H . 
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FIG. 2. Density of states p(E) for the 2d-Andcrson model in the 
presence of resonances with resonance energy Er = 0, coupling 
constant v = 0.2 and various concentrations. The thick dashed 
line corresponds to the Anderson model without resonances. In all 
cases the disorder strength is = 4. 

To verify the validity of the approximations used in the 
derivation of Eq. ^ and to go beyond the weak coupling 
limit, we have performed numerical simulations of the 
2d-Anderson model with m randomly distributed reso- 
nant states. Our main concern in this study is to investi- 
gate how the DOS and the localization length behave as 
a function of energy for various resonance concentration 
/ = m/(N + m) and coupling strengths v. To evaluate 
these quantities, we used the iterative procedure devel- 
oped in jl^. This method is particularly suitable for 
the calculation of various quantities for a macroscopi- 
cally large system described by a microscopic Hamilto- 
nian. The system we considered is a cylinder of width 
M in y direction and length L in x direction. For L 
large compared to M the sample is essentially a quasi- 
one-dimensional system. For the above calculations, we 
have used samples of length L — 10^ and width up to 
M = 100. 

We first examine the behavior of p{E) as a function of 
the coupling constant v and the resonance concentration 
/. Our numerical results for low resonance concentra- 
tion and various values of the coupling v are reported 
in Fig. 1. As we can see, aX E = the DOS shows 
a resonance enhancement which is very well reproduced 
by the analytical formula (^. As the resonance cou- 
pling V increases we observe a gradual decrease of the 
hight of the Lorentzian peak as well as an increase of its 
width. In the limit w = 1 we recover the Anderson case. 
We also investigated the dependence of the DOS on the 
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FIG. 3. Localization length Im{E) in the presence (o) and absence 
(*) of resonances calculated numerically for the 2d- Anderson model 
as a function of energy. The width of the strip is M = 20, and the 
disorder strength is W = 4. The resonance energy is Er = 0, while 
the resonance concentration and coupling are / = 2.4 and v = 0.2 
respectively. It is clear that Im{E) is not affected by the presence 
of resonances. 

Next, we calculate the localization length loa{W,E), 
needed for the determination of diffusion constant D{E). 
Using the same numerical approach pO| we define a lo- 
calization length Im{W,E) depending on M, W and E. 
The asymptotic localization length Iqq{W, E) correspond- 
ing to the infinite system, is then defined in the limit of 
M — > oo i.e. 



1^{W,E)=Im^oo{W,E) 



(8) 



The numerical results, performed for different values of 
resonance concentration, show that l]\i{E) can be consid- 
ered approximately independent from the concentration 
/ and from the resonance coupling v (as long as the con- 
centration is not extremely large). Our results remain 
practically the same for all values of M we have con- 
sidered in our calculations (up to M — 100) and /. In 
Fig. 3 we present the behavior of the localization length 
Im{E) , for one particular value of /, vs. the correspond- 
ing one of the standard Anderson model. This behavior 
of the localization length is totally different from the one 
found in the one-dimensional case where loo (E) decreases 
drastically everywhere except close to the resonance Er 
where it is approximately equal to the corresponding one 



3 



obtained for the standard Anderson model 0. The un- 
derUne physical reason is that the electron in the two- 
dimensional geometry can bypass the resonances follow- 
ing a different path. We expect that the same behavior 
will appear also in three-dimensions. In contrast, in the 
one-dimensional case the electron has to pass through the 
resonance and thus to suffer resonance tunneling. Thus, 
the one dimensional structure acts as a filter allowing 
only electrons with energy Er to transmit. 

We now turn our attention to the calculation of the 
diffusion constant which will be the linking observable 
with the experimental results In the quasi-one di- 
mensional multichannel case the diffusive constant D(E) 
is related with loa{E) and p{E) through the relation |lj] 

loo{E) 



D{E) 



(9) 



Then, it becomes clear from the previous analysis that 
D{E) will show a dip near the resonance energy E,. (see 
also ) , due to the resonance enhancement of the den- 
sity of states p{E). However, in order to compare the 
behavior of D{E) with the experimental results [|| we 
have in addition, to take into account the fact that the 
strength of the resonance coupling v depend on the con- 
centration of resonances / as i; /^/^ 1^. Fig. 4 shows 
D{E) vs. energy for different concentrations. In agree- 
ment with the experimental results Q, we find that an 
increase of / causes a decrease of D far away from the 
resonance while the dip near the resonance energy E^ 
widens and eventually disappears for t; '--^ 1. This result 
identifies the resonant behavior of the DOS as the main 
explanation for the dip in the diffusion constant. More- 
over, it establishes the equivalence of transport proper- 
ties in the presence of resonances, between classical and 
quantum waves. 
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FIG. 4. Normalized diffusion constant D{E^Do{E) calculated 
numerically for the 2D-Anderson model as a function of energy. 
The different curves correspond to different concentrations / of 
the resonances (Do is the diffusion constant in the absence of 
resonances). 



In conclusion, we have studied, both analytically 
and numerically, the effect of resonances on the trans- 
port properties of electrons in two-dimensional tight- 
binding Anderson model. At low concentration and weak 
coupling the DOS is affected strongly by exhibiting a 
Lorentzian peak at resonance energy. As the coupling 
strength increases, the width F of the Lorentzian in- 
creases and the resonance structure in DOS is smeared 
out. Contrary, the localization length is unaffected from 
the resonances. An application of the formula, loo{E) = 



A'Kp{E)D{E) 1 11 1, identifies the behavior of p{E) as the 
explanation for the observed experimental results of the 
diffusion constant. This establishes the analogy between 
classical and quantum waves. Our results are similar to 
those obtained within the ERMT |§. 
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